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Introduction 

In this paper, we give some results concerning the dynamics of double Boolean 
automata circuits (dbac's for short), namely, networks associated to interaction 
graphs composed of two side-circuits that share a node. More precisely, a double 
circuit of left-size £ G N and of right-size r £ N is a graph that we denote by D^ r . 
It has n = I + r — 1 nodes. Nodes that are numbered from to £ — 1 belong to the 
left-circuit and the others plus the node (that belongs to both side-circuits) belong 
to the right-circuit. Node is the only node with in- and out-degree 2. All other 
nodes have in- and out-degree 1. 
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Figure 1: Double circuit D^ r . 



A DBAC D£ tr = (Bg r ,F) is a network denned by (i) its interaction graph, a 
double-circuit D^ r , and (ii) a global transition function F that updates the states of 
all nodes in parallel and that is defined as follows by the local transition functions fi 
of nodes i < n: 

Vx G {0,l} n , 

'F{x)i = fi{x i - 1 ), Vi^{0,£}, 

< F{x)i = fi(x ), (1) 
< F(x) = fo{xe-i,xe+r-2) = fo(xe-i) * fo(x£ +r - 2 ) where * G {A, V}. 

All local transition functions are supposed to be non-constant. Thus, Vi < n, /j, /q, 
/(f G {id, roeg} where id : a h-> a and neg : a i— )• -ia, Va G {0, 1}. As a consequence 
all local transition functions are locally monotone. All arcs (i,j) entering a node j 
whose local transition function is id (resp. neg), with respect to node i, are labeled 
with a + sign (resp. a — sign) and called positive arcs (resp. negative arcs). A side- 
circuit with an even number of negative arcs (resp. odd number of negative arcs) is 
called a positive (side-) circuit (resp. a negative (side-) circuit). 

Given a configuration x = (xq, . . . ,x n -i) G {0, l} n of a DBAC D£ r , we use the 
following notation: 

x L = (x ,. ■ ■ ,x e -i) and x R = (x , x e , . . . , sc n _i). 

A configuration x(t) = (x (t), . . . , x n -i(t)) G {0, l} n such that VA; G N, F k - p (x(t)) 
= x(t+k-p) = x(t) is said to have period p. If x(t) has period p and does not also have 
period d < p, then x(t) is said to have exact period p. An attractor of period p£N, 
or p-attractor, is the set of configurations belonging to the orbit of a configuration 
that has p as exact period. Attractors of period 1 are called fixed points. The graph 
whose nodes are the configurations x G {0, l} n of a network and whose arcs represent 
the transitions (x(t),x(t+ 1) = F(x(t)) is called the transition graph of the network. 

In [T], the authors showed the following results: 

Proposition 1 1. The transition graphs of two DBACs with same side-signs and 
side-sizes are isomorphic, whatever the definition of /q (i.e., whether ★ = V or 
* = A in the definition |ip of F above). 

2. Attractor periods of a DBAC divide the sizes of the positive side-circuits if there 
are some and do not divide the sizes of the negative side-circuits if there are 
some. 
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3. If both side- circuits of a DBAC Di^ r have the same sign, then, attractor periods 
divide the sum N = £ + r. 

4- If both side-circuits of a DBAC Dm have the same sign and size, then Di^ 
behaves as an isolated circuit of that size and sign (i.e., the sub-transition graph 
generated by the periodic configurations of Dj>g is isomorphic to the transition 
graph of an isolated circuit with the same sign and size). 

5. A DBAC has as many fixed points as it has positive side- circuits. 

6. If both side-circuits of a DBAC Dg r are positive andp divides t and r, then num- 
ber of attractors of period p is given by A+ (sequence A1037 of the OEIS f^), 
namely, the number of attractors of period p of an isolated positive circuit of 
size a multiple of p. 

As a result of the first two points of Proposition [TJ we may focus on canonical 
instances of DBACs. Thus, from now on, we will suppose that ★ = V and \/i ^ 
0, fi = id. If the left-circuit is positive (resp. negative), we will suppose that the arc 
(£ — 1,0) is positive (resp. negative) and f$ = id (resp. f$ = neg) and similarly for 
the right-circuit. Thus, the only possible negative arcs on the DBACs we will study 
will be the arcs {I — 1,0) and (n — 1,0). 

The last point of Proposition [T] yields a description of the dynamics of a doubly 
positive DBAC (in terms of combinatorics only but [Tj also gives a characterisation 
the configurations of period p, Vp G N for all types of DBACs). Thus, here, we will 
focus on the cases where the DBACs have at least one negative side-circuit. 

For a negative-positive or a negative- negative DBAC Di^ r , we will use the following 
notations and results. The function fi is the Mobiiis function. It appears in the 
expressions below because of the Mobiiis inversion formula [2j. 

• The number of configurations of period p is written: 

C p (£,r). 

• The number of configurations of period exactly p is written and given by: 

c;(£,r) = £>(^)-C p (£,r). 
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• The number of attractors of period p (where p is an attractor period) is written 
and given by: 

W) = = -■!>(-) -Cp&r). 

q\p 

In particular, Ct(£,r) = C$(£,r) = 1. 

• The total number of attractors is written and given by: 



E 

p attractor period 



A p (£,r) 



For the sake of clarity, we will write A p (£,r) = A^(£,r) (resp. A p (£,r) = A^{£,r)) 
and T(£,r) = T^(£,r) (resp. T(£,r) = T = (£,r)) when L>^ r will be a negative-positive 
DBAC (resp. a negative-negative dbac). 

1 Positive-Negative 

We first concentrate on DBACs whose left-circuit is negative and whose right- 
circuit is positive. From Proposition [T] we know that all possible attractor periods 




Figure 2: Interaction graph of a negative-positive DBAC. All arcs are positive (resp. all 
local transition functions are equal to id) except for the arc {£ — 1,0) (resp. except for the 
local transition function Jq — neg). 

of such DBACs divide r and do not divide £. We also know that these networks have 
exactly one fixed point. In the sequel, we focus on attractor periods p > 1. 
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Table 1 : Total number of attractors of a negative-positive DBAC Dt^ r (obtained by computer 
simulations). Each colour corresponds to a value of gcd(£, r). The line (resp. the column 
TJ) gives the total number of attractors of an isolated positive (resp. negative) circuit of 
size r (resp. £). 
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Characterisation of configurations of period p 

Let p be a divisor of r = p ■ q that does not divide £ = k ■ p + d, d = £ mod p. 
We write A p = gcd(p,£) = gcd(p,d). Let x = x(t) be a configuration of period p: 
Vm G N, x(t + m-p) = x(t). Note that to describe the dynamics of Di^ r , it suffices to 
describe the behaviour of node (see [1]). The configuration x satisfies the following: 

xo(t) = x (t + k-p) 

= -<X£-i(t + k ■ p — 1) V x n -i(t + k ■ p — 1) 
= ^x e _ k . p (t) V x (t + (k - q) ■ p) 
= ->Xd(t) V xo(t) 

Thus, if Xo(t) = 0, then Xd(t) = 1 (and also, because xo(t) = ->xo(t — d) V 
xo(t), if xo(t) = 0, then xo(t — d) = 1). From this, we may derive the following 
characterisation: 

Proposition 2 Let p G N be divisor of r = p- q that does not divide i. A configura- 
tion x = x(t) has period p if and only if there exists a circular word w £ {0, 1} P of 
size p that does not contain the sub-sequence OuO, u G {0, l} ' -1 and that satisfies: 

x L = w k w[0...d- 1] and x R = w q 

where w[0 . . .m] = wq . . . w m . More precisely, the word w satisfies: 

Vi < p, Wi = xo(t + i). 

Consequently, the property for a configuration to be of period p depends only on 
d = i mod p and on p (and not on £ nor on r). 

As a result of the characterisation above, we may focus on negative-positive 
DB ACs Di )T such that t < r ( i.e., £ = £ mod r) and to count the number of attractors 
of period p we may focus on DBACs Dg T such that £ < p (i.e., £ = d = £ mod p). In 
other words, from Proposition [2] 

(£, r) = (£ mod r, r) = Ap (£ mod p, r) . 

Combinatorics 

From Proposition [2] each configuration x of period p is associated with a circular 
word w G {0, 1} P that does not contain the sub-sequence OuO, u G {O,!} ^ 1 , d = 
£ mod p. It is easy to check that this word w can be written as an interlock of a 
certain number iV of circular words w/ 1 ),?// 2 ), . . . , w( N ' of size T = p/N that do not 
contain the sub-sequence 00 (see figure [3|. More precisely, the size of a word 
satisfies T -d = K -p for a certain K such that T and K are minimal. In other words, 
T ■ d = lcm[d,p) = gJiu p ^ and thus T = Consequently, iV = A p and we obtain 
the following lemma which explains why, in each column of Table [TJ every cell of a 
same colour contains the same number. 




Figure 3: The circular word w — wq...w p -i = Xo(t) . . . Xo(t + p — 1) mentioned in 
Proposition [2] that characterises a configuration x{t) of period |tu| = p. In this example, 
p = 15, d — I mod p — 6 so that w is made of an interlock of A p = gcd(d,p) = 3 words 
10 W = Wq . . . 104 of size p/Ap = 5. 

Lemma 1 Lei p be a divisor of r = p ■ q that does not divide t. The number of 
configurations of period p, C p (£,r), depends only on A p = gcd(£,p) and on p. Thus, 
we write: 

C p (£,r) = C pAp . 

The number of circular words of size n that do not contain the sub-sequence 00 
is counted by the Lucas sequence (sequence A204 of the OEIS [3]): 

'L(l) = l 
< L(2) = 2 

L{n) = L(n - 1) + L(n - 2) = 0" + T = + 

where <f> = ~ 1.61803399 is the golden ratio and 4> = 1 - <f> = ~ 

—0.61803399. Among the properties of (j) that will be useful to us in the sequel are 
the following: 

(f 2 = 1 + <p and (p = . 

4> 

Thus, to build a circular word w £ {0, 1} P without the sub-sequence OuO, u 6 
{0, l}^ 1 , one needs to chose A p among L(-f-) words w^' of size without the 
sub-sequence 00. As a result holds Proposition [3] below: 
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Proposition 3 The number of configurations of period p is given by: 

Consequently, the number of attractors of period p is given by: 

A p ± (^r)=A p , Ap = i.^^).L(^)^, 
P q \ P q 

where A p = gcd(£,p). 

Number of configurations of period p 

Let us develop the expression for C p (£,r) = C p ^ p : 

C P A P = £(£) Ap 

= (-0r p -E fc <A P (t)-^'^-(- 1 )^ 

= 4? 

= 4? ■ ((-(A 2 )^ + l) AP 

= (-l)p- ((-1)^ • (0 2 )^ + 1) A P. 

If p is odd, cannot be even. Thus, there are three cases only: 

1. p and are odd. Thus, because A p is necessarily also odd: 

c pAp = -MM-(<A 2 )* + i) Ap = HM(<A 2 )*-i) Ap . (2) 

2. p is even and is odd. Thus, because A p is necessarily even: 

c pAv = HM-W> 2 )* + i) Ap = H p -((0 2 )£-i) A '. (3) 

3. p and are both even. Thus: 

C pAp = |^-((0 2 )4 + l)A P . (4 ) 
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To sum up, we give below Proposition [4] whose last part can be derived from 
the relation between the Euler totient and the Mobiiis function /u(-), ip(n) = 
Ylm\n( n / m ) ' M 771 ); an d from the following equations where A q = gcd(q,£): 

p\r p\r q\p ' p\r q\p 

q\r ' q\r 



Proposition 4 Let A p = gcd(£,p) where p is a divisor of r that does not divide £. 
Then, the number of configurations of period p is given by: 



C pA p 



■ P ■ 



p 

h2 \ A„ 



l) Ap if ls odd ' 



h2\ A 



«• +1) 



is even. 



In particular, C r ^ c ^^ r ) counts the total number of periodic configurations of the 
network. The number of attractors of period p and the total number of attractors are 
respectively given by: 



-•yV-J-C^A,, arid T±(^r) 
p ' q 

i\p 



p\r, -*(p\£) 



Upper bounds 

From equations Q, (j3j) and Q, one can derive that C pAp and thus (£, r) = 
^pA are max i ma l when A p is minimal (i.e., A p = 1), if p is odd and if p is even, 
on the contrary, C pAp and h p (£, r) are maximal when A p is maximal (i.e., A p = |). 
Thus, we have: 



< |^|P.((^) + i)f = Cp P_ 



(3+3-. 



4>p 

( 3+3-0 

p 



32, 
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In addition, < i • F p (y/3) where F p (a) = Sd|p /■*(«?) ' Now, it can be shown 

that: 

Va > 0, Vp / 2, a v ~ x < F p (a) = ■ ^ < a p . 

dip 

Consequently, 

Vp + 2, A p , Ap < 1 • F P (V3) < 2 • {^f ■ \ ■ F p {2) = 2 • A* • A+ (5) 

where A+ refers to the number of attractors of period p of an isolated positive circuit 
whose size is a multiple of p [I]. Then, from (pi), we may derive the following result: 



Proposition 5 For all p ^ 1, the number of attractors of period p satisfies: 

4(£,r)=A pAp < -A+ 
and for p = 2: A2 5 a 2 = 1 = A^\ 

Let us denote by the number of attractors of a positive circuit of size n. From 
Proposition [5j for all r, the total number of attractors of a negative-positive DBAC 
satisfies: 




However, computer simulations (c/ Table ll last line) show that this bound on 
is too large. We leave open the problem offinding a better bound. 



2 Negative-Negative 

We now concentrate on doubly negative DBACs. The canonical DBAC we will use 
in the discussion below is defined in Figure [4] 

Let p £ N be a possible attractor period of Di jr {p divides N = £ + r but divides 
neither £ nor r). Without loss of generality, suppose £ mod p > r mod p = d. 
Because p divides £ + r, it holds that £ mod p = p — d. Then, for any configuration 
x = x(t) £ {0, l} n of period p, we have the following: 

x (t) = ->X£-x(t- 1) y-ix n -i(t— 1) 

= -ix (t - £) V -ix Q (t - r) 

= -ixolt + r) V^x (t + £) 

= -ix (t + d) V-ix (t- d) 

As a consequence, if Xo(t) = 0, then xo(t + d) = xo(t — d) = 1 and if Xa(t) = 
1, then either xo(t + d) =0, or xo(t — d) =0. Thus, the circular word w = 
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Figure 4: Interaction graph of a negative-negative DBAC. All arcs are positive (resp. all 
local transition functions are equal to id) except for the arcs (£ — 1,0) and (n — 1,0) (resp. 
except for the local transition functions ffr = neg and = neg). 

Xo(t) . . . xo(t + p — 1) contains neither the sub-sequence OuO nor the sub-sequence 
lulu'l (u,u' € {0,l} d_1 ). 

Let A = gcd(£,r). As in the previous section, w can be written as an interlock 
of A p = gcd(d,p) = gcd(A,p) words w^' of size p/A p that do not contain the sub- 
sequences 00 and 111. As one may show by induction, the number of such words is 
counted by the Perrin sequence sequence A1608 of the OEIS [3j: 

>(0) = 3, 
( P(2) = 2, 

P(n) = P(n - 2) + P(n - 3) = a n + (3 n + f3 n , 

where a, j3 = ^ ■ (—a + % ■ ^ — 1) and (3 are the three roots of x 3 — x — 1 = 0, and 
a, the only real root of this equation, is called the plastic number [5]. 

Using similar arguments to those used in the previous section, we derive Propo- 
sition [6] below. This proposition explains why, in Table [3j all cells of a same diagonal 
(i.e., when = I + r is kept constant) that have the same colour also contain 
the same number: the number of attractors depends only on JV = I + r and on 
A = gcd(£,r) and not on i nor r. Equations in Proposition [6] exploit, in particular, 
the fact that if p divides I or r then C Pj a p = (because then P(-£-) = P{1) = 0). 

Proposition 6 Let N = i + r and let p € N be a possible attractor period of D^ )T (p 
divides N but divides neither I nor r). Let also A = gcd(l,r) and A p = gcd(A,p). 
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Then, the number of configurations of period p of the doubly negative DBAC -D^ r 
depends only on p and A p . It is given by: 

C P (£,r) = C pAp =P(^-) A v. 

The number of p-attractors and the total number of attractors of a doubly negative 
DBAC Di^ r are respectively given by: 




The expression for T^y A in Proposition [HI above simplifies into the following if 

_ N ■ : ' 



K = ^ is a prime: 



N 

q\A, gcd(q,K)=l 

In particular, ifi^=^=2orET=^=3, then because P(2) = 2 and P(3) = 3, 
the following holds: 



g|f , fl od(ff,2)=l 



1 * — T N 

3 



q|f , sai( g ,3)=l 



From the computer simulations we performed (see Tables [2] and [3]) , we observe 
the following: 

1. Given £, T = (£, r) is maximal when r = I, 

2. Given an integer N which is not a multiple of 3, is maximal when A 
is maximal (in particular, if N is even without being a multiple of 3, then 

T iV,A — T /v s) m 

> 2 

3. Given an integer N which is a multiple of 3, T^- A is maximal when A = y. 
We leave the proofs of these three points as an open problem. 
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Table 2: Total number of attractors of a negative-negative DBAC Dt^ r (obtained by computer 
simulations) . 
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33 


19 


23 


28 


278 


46 


60 


75 


88 




— ^ 


2048 


17 


11 


11 


16 


19 


24 


28 


39 


46 


60 


75 


97 






3856 


18 


11 


17 


17 


23 


28 


6 


46 


60 


729 


96 






7286 


19 


16 


19 


24 


28 


39 


46 


60 


75 


97 






13798 


20 


19 


23 


28 


32 


125 


60 


75 


88 






26216 


21 


24 


28 




46 


60 


96 






— V 


49940 


22 


28 


38 


46 


60 


75 






95326 


23 


39 


46 


60 
66 


75 


97 






182362 


24 


46 


60 


88 




X 


349536 


25 


60 


75 


97 




£ = 4 x r 


671092 


26 


75 


96 






1290556 


27 


97 






2485534 



Table 3: Total number of attractors of a negative-negative DBAC D( r (obtained by computer 
simulations). Each colour corresponds to a value of gcd(£, r). The last column gives the total 
number T7 of attractors of an isolated negative circuit. 
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